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Abstract 

We employ an automated detection algorithm to perform a global study of solar 
prominence characteristics. We process four months of TESIS observations in 
the He ll 304 A line taken close to the solar minimum of 2008-2009 and focus 
mainly on quiescent and quiescent-eruptive prominences. We detect a total of 
389 individual features ranging from 25 x 25 to 150 x 500 Mm^ in size and obtain 
distributions of many their spatial characteristics, such as latitudinal position, 
height, size and shape. To study their dynamics, we classify prominences as 
either stable or eruptive and calculate their average centroid velocities, which 
are found to be rarely exceeding 3 km s“^. Besides, we give rough estimates 
of mass and gravitational energy for every detected prominence and use these 
values to evaluate the total mass and gravitational energy of all simultaneously 
existing prominences (lO^^-lO^^ kg and lO^^-lO^^ erg, respectively). Finally, we 
investigate the form of the gravitational energy spectrum of prominences and 
derive it to be a power-law of index —1.1 ± 0.2. 

Keywords: Prominences, Quiescent; Prominences, Dynamics; Prominences, 
Formation and Evolution 


1. Introduction 


Prominences are one of the most noticeable features of the Sun, which, although 
observed for over a century, are still far from being completely understood. 
For historical reasons, a distinction is made between prominences, observed 
off-limb as luminous formations, and filaments, usually seen in absorption on 
the disk. Physically, these are the same structures consisting of plasma with 
properties similar to those of the chromosphere being nearly 100 times denser 
and cooler than the surrounding corona (Hirayama, 1985 Tandberg-Hanssen, 
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1995 Labrosse et al.^ 2010). Prominences vary largely in their size, morphology, 


dynamics, magnetic structure and the way they are related to the rest of the 
solar activity. As more of these differences were revealed, new classification 
schemes were introduced by a number of authors (Pettit, [1932 Severny and 
Khokhlova, 1953 Zirin, 1966 Leroy, Bommier, and Sahal-Brechot, 1984). In this 
paper, we will follow the most general one, which divides prominences into two 
major classes depending on whether they are part of an active region (active 
prominences), or whether they exist separately over the quiet Sun (quiescent 
ones). It is known that instabilities in the magnetic structure supporting both 
these types of prominences can result in eruption of some part of their material 
high into the corona (van Ballegooijen and Martens, 1989 Zhang, 2013). To such 
prominences we will refer here as eruptive ones. 

Comparative studies have shown that more than 70% of prominence eruptions 
are followed by coronal mass ejections, or CMEs (Gilbert et a/., |2000| Gopal- 


swamy et a/., 2003 Filippov and Koutchmy, 2008| ). They, in turn, travel over 
vast distances from the Sun and have a considerable impact on the heliosphere 
in general, including the Earth’s magnetosphere. At the same time, the exact 
connection between prominence eruptions and GMEs, as well as the mecha¬ 
nisms of formation and eruption of prominences themselves, are still poorly 
understood. Existing models of prominence formation consider a number of 
possible magnetic topologies and mechanisms that can fill them with plasma, 
while observational data are not sufficient to discriminate between these models 
(Galsgaard and Longbottom, 1999 Mackay et a/., 2010). Similarly, theories 
of prominence eruptions developed so far suffer from poorly constrained pre¬ 
emptive conditions of the magnetic field and plasma (Lin, Soon, and Baliunas, 


2003). In this context, statistical studies of prominences seem particularly useful 


as they may provide the necessary conditions and constraints to theories and 
establish reliable correlations of prominence evolution with other manifestations 
of solar activity. 

For decades, prominences have been observed from the Earth in a num¬ 
ber spectral lines, primarily as well as in radio emission (Labrosse and 
Gouttebroze, 2001 Heinzel, 2007). A comprehensive study of these phenomena, 
however, became only possible with the advent of space-based observations. 
Nowadays, prominences can be observed from space by a number of instruments 
and in a much wider spectral range. Above all, these observations are free of 
day-night interruptions and unfavourable atmospheric effects. Most of them are 
performed in the EUV spectral range, where prominences are best seen in the 
He II 304 A line. Despite the vast amounts of data produced by these instru¬ 
ments, most of the present-day studies of prominences limit themselves to only 
examining individual, and usually, the largest events. In this way, however, one 
risks missing smaller prominences, which are of great interest still, and obtaining 
an incomplete picture of these phenomena. Therefore, to make use of all the 
data available, and to study the entire set of prominences, special detection and 
tracking algorithms need to be applied. 

A continuous set of observations produced by the Extreme ultraviolet Imag¬ 
ing Teleseope (FIT) of the Solar and Heliospherie Observatory (SOHO), the 
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Extreme Ultraviolet Imager (EUVI) of the Solar TErrestrial RElations Obser¬ 
vatory (STEREO) and the Atmospherie Imaging Assembly (AIA) of the Solar 
Dynamies Observatory (SDO) since 1996 now covers a vast period of almost 
two solar cycles. Nevertheless, only a few attempts have been made so far 
to use these data for large statistical studies of prominences. Two algorithms 
have been developed for EIT data (Eoullon and Verwichte, 2006 Labrosse, 


Dalla, and Marshall, 2010) and another one for STEREO-A and -B observations 
(Wang et al.^ 2010). Eurthermore, only the last attempt produced a prominence 


catalogue that is now publicly available on-line. Continuing this line of study, 
we have developed a new algorithm capable of precisely locating prominences 
using monochromatic observations in the He ll 304 A line only. We employ this 
algorithm to study a period of four months of observations, which corresponds 
roughly to the end of the prolonged solar minimum of 2008-2009. Hence, we 
focus mainly on quiescent and quiescent-eruptive prominences and leave active 
prominences for future studies. 

The organization of this paper is as follows: in Section we present the 
observational data used and in Section we give a general description of the 
detection algorithm employed. We then review the major results of its operation 
and perform statistical analysis of the obtained prominence parameters in Sec¬ 
tion]^ Einally, we discuss the system’s performance, as well as advantages and 
future prospects of automated detection of prominences in Section 


2. Observations 


Eor this study, we use observations from one of the two Eull-disk EUV Tele- 
seopes (EET) of the TESIS solar observatory, namely, the EET-304/171 tele¬ 
scope (Kuzin et a/., 2009). This instrument was designed at the Lebedev Institute 
of the Russian Academy of Sciences and launched on board the Complex ORbital 
Observations Near-Earth of Aetivity of the Sun-Photon (CORONAS-Photon) 
satellite on 30 January 2009. To our advantage, CORONAS-Photon’s relatively 
short operation time (until mid January 2010) coincided with the end of an 
extraordinarily long solar minimum between Solar Cycles 23 and 24. This gave us 
an opportunity to unrestrictively observe small prominences, barely seen during 
the rest of the cycle against larger manifestations of solar activity. Here, we 
analyse the most complete portion of this dataset from August to November 
2009, which is publicly available on-line at http://tesis.lebedev.ru, 

With a field of view of 1°, the telescope was able to image the full solar disk 
and the corona up to the distance of ^ 0.95 solar radii above the limb (see an 
example of a typical observation in Eigurej^. It produced 2048x2048 images 
with a pixel size of ^ 1.7 arcsec, or roughly 1.25 Mm. Half of the images, however, 
have been reduced to 1024x1024 due to telemetry limitations. Eor the regular 
observations in the 304 A channel the cadence was set to 4.0 min, but as the 
observations were non-continuous because of the partly Earth-shaded orbit of 
the satellite, the resulting mean cadence amounted to ^ 7.2 min. To make sure 
that the CCD showed no apparent signs of degradation during the observation 
period, which could affect our investigations, we calculate the average response 


SOLA: promstat.tex; 2 July 2015; 1:00; p. 3 













I.P. Loboda and S.A. Bogachev 



Figure 1. Example of typical TESIS observation in He ll 304 A line taken on 26 September 
2009 at 15:44:46 UT (intensities are in logarithmic scale). Note that solar North is not up in 
the original images. 


of the CCD, shown in Figure It clearly shows that possible degradation of 
the CCD is negligible against the larger variations in solar irradiation. 

Since the entrance and detector filters of the telescope were relatively broad, 
dedicated mainly to block the strong white-light radiation from the Sun, the 
bandpass in the 304 A channel of TESIS was predominantly defined by the 
spectral reflectivity of its focusing mirror, which is shown in Figure (Kuzin 
et a/., |2011|). This figure also indicates the main contributors to this relatively 


broad-band channel: the closely spaced Si xi line at 303.3 A and He ll line 
at 303.8 A. Contribution of the Si xi line has been previously investigated by 
several authors, who have found it to be around 3-4% for the on-disk quiet 
Sun regions and more than 90% for the coronal emission (Cushman and Reuse, 


1978 Thompson and Brekke, 2000). Since plasma conditions in prominences are 


similar to that of the quiet Sun, helium emission is typically dominant in them. 
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(a) 



(b) 



Figure 2. |(a)| Averaged signal from the CCD throughout the observation period; we have 
made the correction for varying Sun-Earth distance. |(b)| Spectral reflectivity of the focusing 
mirror in the 304 A channel of TESIS and the corresponding EUV lines. 


3. Image Processing 

A total of 18812 valid images of the Sun in the He ll 304 A line were ob¬ 
tained during the specified period. As we mentioned above, processing such an 
amount of data, including detection of prominences and measurement of their 
parameters, is not feasible by hand. For this purpose, we have developed an 
automatic method, capable of locating the prominences on individual 304 A 
images, tracking these features throughout their lifespan, and determining their 
most important spatial and dynamic characteristics. The operational principle 
of this method is described below as a sequence of several processing steps. 


3.1. Preprocessing and Background Removal 


In the preliminary stage, the data is prepared for subsequent operations: cor¬ 
rupted images are filtered out and for the remainder corrections for minor 
defects are made (e.^., hot pixels and cosmic ray hits). Apart from that, off- 
limb regions of the images are transformed to polar coordinates and stored as 
rectangular arrays in order to simplify further image processing (for an example, 
see Figure 4a). 


Essentially, the next step would be to detect prominences on these pre- 
processed images. Unfortunately, as one can see from the Figureprominences 
are observed against a strong background, which can be attributed either to the 
thermal noise of the CCD and extra signal due to the scattered radiation, or to 
the radiation of the undisturbed corona itself. As Labrosse, Dalla, and Marshall 


(2010) rightly note, to know the background is critical to the whole process. 
Analysis of individual images shows, that this background does not change 
much with time, but is spatially highly inhomogeneous: it rapidly decreases with 
height and varies with latitude, having significant dips around the poles. To our 
advantage again, because the observation period is close to the solar minimum, 
the positions of the polar coronal holes remain stable in spite of the solar rotation. 
This allows us to determine the background following the statistical approach 
described below. 
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(a) 


(b) 



Figure 3. | (a) [ Radial profiles of the quiet Sun (solid grey line), prominence (solid black line), 
statistically obtained background (dotted line), anal ytical fit to that background (dashed line) 
and its exponential part (dashed-dotted line). [(1^ Intensity distribution at a fixed point of 
corona (grey line) with an assumed Gaussian contribution from the quiet Sun (black line). 


First, we use the whole dataset to build the intensity distribution of every indi¬ 
vidual pixel of the previously polar-transformed and aligned images, an example 
of which is shown in Figure [31^ We assume that in the quiet Sun, the observed 
intensity coming from a fixed point in the corona follows a normal distribution. 
This quiet Sun distribution is supplemented with stronger signals coming from 
a small number of active features, which occasionally emerge in the corona. 
Therefore, we approximate the left slope of the obtained distribution with the 
Gaussian function and take the position of its maximum as the estimated value 
of background in that pixel. 

In order to smooth out the noise and small-scale irregularities, we then look 
for an analytical function to fit the background we have thus obtained. Having 
studied its radial and latitudinal profiles on semi-log and log-log plots, we have 
finally arrived to the following model: 

h) = , (1) 

where h is the height above the limb, is the heliographic latitude and a{ip)^ 
b{ip), c{ip) and d{ip) are generalized logistic functions of a similar form, e.g. 




a2 — CLl 




( 2 ) 


In other words, a{(p) is equal to ai at the equator and middle latitudes and to 
a 2 in the vicinity of poles, with a smooth transition from ai to a 2 centred at 
the latitude (/^o, the steepness of this transition determined by 5ip. Parameters 
ifo and 5if are taken the same for all four functions. 

Although the double exponent in Equation 0 may seem strange at a first 
glance, it proved to be the best approximation to the specific radial profiles of 
the background obtained for this dataset; while B(h) has an almost exponential 
dependence at higher altitudes, it becomes much steeper close to the limb (Fig- 
3a). Moreover, below a height of ^ 15 Mm the background is mostly due to 


ure 


the strong radiation from the EUV spicules and prominences are in general not 
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Figure 4. Processing steps of the image in Fignre | (a) | polar-tr ansfo rmed image with minor 
defects corrected, latitudes counted clockwise from the North pole, [(^background subtracted, 
1(c) I normalized image, [(d)] kernels used for the subsequent region-growing procedures, |(e)| final 
selections. Intensities in panels a to c are in logarithmic scale. 


visible there. Therefore, we do not attempt to fit this part with an analytical 
function. 

Subsequently, we follow a technique similar to that described in Morgan, 
Habbal, and Woo (2006) and use this modelled background B to normalize 
every image in the dataset according to the equation: 


I^S 




B^6 


I - 


(3) 


where / = I{(p,h) is a given polar-transformed image and N = N{ip^h) is the 
same image after normalization (see an example in Figure [^. We introduce here 
a small parameter 6, first, to avoid division by zero, and second, to suppress the 
noise at higher altitudes, where B is close to zero. Thus, it must be substantially 
greater than the average noise level in the periphery of the image 5 CCD 
counts for this instrument) and substantially smaller than typical signal values 
from prominences 100 and more CCD counts). For the dataset employed, we 
have set 6 equal to 25. 


3.2. Prominence Detection in the Individual Images 

To avoid restrictions imposed by the necessity of simultaneous multi-wavelength 
observations, our first concern was to be able to detect prominences using 304 A 
observations only. The problem is, however, that they are not the only bright 
structures seen off-limb in this usually broad-band channel. High coronal tem¬ 
peratures above the active regions (ARs) can stimulate emission in the Si xi line 
and cause a considerable increase in brightness above the limb. Consequently, 
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(a) 


(b) 




Intensity threshold Tj 

Figure 5. Sele cted area Ai as it grows with the decrease of the cut-off threshold Ti in case 
of a prominence | (a) | and an active region [(b)] 


additional procedures are required to discriminate between true prominences 
and these AR-imposed effects. A number of solutions to this problem have 
been previously proposed by several authors. Foullon and Verwichte (2006), for 
example, implemented a histogram segmentation to set an intensity threshold 
for an individual image and then employed multi-wavelength observations to 
discard active regions. Wang et al. (2010) first selected potential prominences 
by a preset threshold and then performed a regression analysis based on the 
shape of the selected regions. A disadvantage of this method is, however, that 
wrong thresholding may lead to misclassification and therefore to no prominence 
detection. Finally, an interesting method was proposed by Labrosse, Dalla, and 
Marshall (2010); these authors used moments to find radial profiles character¬ 
istic of prominences only and then reconstructed their shape by morphological 
opening procedures. 

We have developed a different technique, taking advantage of the fact that, 
in the He ll 304 A line, quiescent prominences tend to have sharp boundaries 
(at least in the dataset concerned), and on the contrary, active regions reveal 
themselves off-limb as diffuse brightenings. In our study, we concentrate mainly 
on the measurement of the physical parameters of prominences for the later 
statistical treatment. Thus, our primary concern was to achieve high precision 
in prominence detection in order to derive their properties most reliably. In 
order to meet this requirement, our algorithm individually determines an optimal 
threshold for each prominence on the image by the following procedure. 

First, kernels of future prominences are found to isolate individual off-limb 
features. These are the areas where the normalized signal N is greater than a 
fixed threshold Tkem = 3.0. Each kernel is then treated separately and expanded 
by means of the region-growing procedure. This operation is performed multiple 
times with a varying cut-off threshold T^, which is gradually decreased from 
Tmax = 7.0 to Tmin = 1-1 ^t an increment of ST = —0.02. Concurrently, the 
area of the resulting selection Ai is measured as a function of T^. Naturally, as Ti 
decreases, Ai grows monotonically; a typical view of this dependence is shown 
in FigureOne can notice that Ti can be greater than Tkem- In this case, the 
region-growing procedure is reduced to finding all the fitted points inside the 


SOLA: promstat.tex; 2 July 2015; 1:00; p. 8 









Quiescent and Eruptive Prominences at Solar Minimum 



Figure 6. Resulting selection at several stages of the iterative region-growing procedure. 


kernel. In addition, if in the region-growing process two kernels merge together, 
their measured dependencies Ai(Ti) are summed. 

Initially, when the main body of a prominence is being selected, Ai grows quite 
rapidly (Figures [^pj|d|). Later, when the sharp boundaries of the prominence are 
reached, the increase of Ai slows down (Figures [^p[|h|); this is easily identifiable 
in Figure [5^ as some sort of a saturation plateau. Finally, when becomes suf¬ 
ficiently low such that regions surrounding the prominence begin to be selected, 
Ai bursts to a rapid growth (Figures [^|[|I|). The value of Ti preceding this burst 
is then chosen as an optimal intensity threshold for the given prominence. 

The main advantage of this approach is the high dynamic range of detection. 
This is well illustrated by Figure at the late stages of eruption the faint 
prominence material is detected just as well as the bright body of the prominence 
before the eruption has started. Another important feature of this technique is 
that it allows the direct discrimination between prominences and active regions. 
As the latter have no sharp boundaries, the dependence Ai{Ti) is smoother, with 
no characteristic knee as is the case of a prominence (Figure [51^. Therefore, no 
optimal threshold is determined and they are not detected as prominences. This, 
of course, inhibits the detection of active prominences, that are often inseparable 
from the active regions, but as we state above, they are not the subject of the 
present study. 

We should note, however, that this technique is still biased by the arbitrary 
choice of Tkem- Fortunately, this parameter is chosen such that it would fall into 
the saturation plateau in Figure [5a| which precedes the burst of the selected area. 
Since the plateau is broad enough, there is no difficulty to choose T^em^ which 
will fall into this interval for most of the prominences. Therefore, we believe that 
small variations of T^em a-re not relevant for the algorithm performance. When 
testing the algorithm, we have eventually chosen Tkem = 3.0 since it produced 
better results for a small portion of the dataset that we have studied manually. 
As for the computational efficiency of the algorithm, the most time-consuming 
operation is the iterative region-growing procedure described above. However, 
since each successive iteration does not expand the initial seed, but continues to 
expand the result of the previous iteration, the whole process can be interpreted 
as a single process, with a number of intermediate stops at certain thresholds. 
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Figure 7. The high dynamic range of detection exemplified by a large eruptive event observed 
on north-eastern limb on 26 September 2009. Result of the detection outlined with a bright 
curve; intensities are in logarithmic scale. 


Therefore, the actual processing time is much smaller than that required to 
preform the same number of independent region-growing procedures. In our case, 
the processing of each image took typically from 2 to 3 seconds on a 3.1 GHz 
eight-core CPU, depending on the number of kernels initially selected. 

3.3. Parameter Extraction, Event Tracking, and Catalogue 
Compilation 

Once all prominences are duly located on the image, for each of them, helio¬ 
graphic latitude and altitude of the highest point are determined, and visible (z.e. 
projected) area, integral brightness, and position of the centroid are calculated. 
Furthermore, a set of models described below in Section is used to estimate 
its volume, mass and, therefore, its gravitational energy. Additionally, in order 
to find the erupted parts of prominences, selections that are detached from the 
disk and are entirely located above a certain height of 30 Mm are marked. 

Eventually, the last step of our algorithm is to track the same prominences as 
they appear on a series of consecutive images. We notice that most of the promi¬ 
nences are confined to the same latitudes during their entire lifetime. Therefore, 
two prominences can be considered as the same event if they overlap on two or 
more consecutive images. They can be absent, however, on no more than two 
images in a row and must be observed on a minimum of three images. Having 
tracked a prominence with this approach, the algorithm registers the moment 
of its appearance and disappearance, calculates its lifetime and the number of 
images in which it was detected. For all of the characteristics of a prominence 
determined previously on the individual images, the minimum, the maximum, 
the average, and the standard deviation are found. After that, tangential and 
radial velocities of the centroid of a prominence are calculated based on the 
change of its position, as well as the radial velocity of the highest point of the 
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prominence. Finally, if the tracked prominence had at any time regions that 
were marked as erupted, the whole event is classified as eruptive prominence. 
All of these data are stored in a table, which is the final result of the application 
of the algorithm. In the next section, we briefly describe these results for the 
observation period and perform a statistical analysis of the obtained properties 
of the prominences. 

To evaluate the accuracy of the detection, we have manually investigated a 
relatively short period of ten days from 10 to 20 September 2009. During this 
period, a total of 46 prominences was detected. We found observable promi¬ 
nences in the images for 42 of these detections, although some of them were 
barely distinguishable by eye. For the remaining four, we could find no features 
resembling an off-limb prominence. Most likely, these false detections resulted 
from the increased coronal radiation due to the presence of an internetwork 
bright point close to the limb, or from a superposition of short-lived and small- 
scale phenomena, like surges and polar macrospicules. At the same time, we have 
found no large prominences that were not detected by the algorithm; it is still 
possible, however, that we could have missed several rather small ones. Besides, 
all 42 prominences were correctly classified by the algorithm as either stable or 
eruptive (although, for smaller prominences it proved to be harder to draw this 
distinction). The algorithm was sometimes unable to detect faint remainders or 
precursors of the prominences, therefore the derived lifetimes are most likely to 
be slightly underestimated. 


4. Results 

4.1. Dynamics and Morphology 

For the period of four months, we have detected a total of 389 quiescent promi¬ 
nences, 41 of them were classified as eruptive. Detected prominences vary largely 
in size, the biggest is 500 Mm long and 150 Mm high and the smallest 25 Mm both 
in length and height (Figures |8a|b| ). Thus, the characteristic size of a prominence 
L = ^/LH (where H is the prominence height and L is the prominence length 
along the limb) is less than 70 Mm for 50 % of the prominences. Most of the 
detected prominences are stretched along the limb; for only 10 % the ratio L/H 
was found to be less than 1 (Figure [^. As for the latitudinal distribution of 
prominences, more than 95 % are concentrated at latitudes below ±60°, the 
maximal probability density falls at medium latitudes around ±40° (Figure [8d]). 

Most of the detected prominences proved to be stable, showing practically no 
apparent bulk motion. This is well demonstrated by Figure for more than 
92 % of the prominences the centroid velocity, averaged over the prominence 
lifetime, does not exceed 3 km s“^. This value, however, does not reveal the 
inner plasma motions in prominences. One of the surprising results is that around 
75 % of the prominences were detected off-limb for less than 24 hours (Figure [^. 
Manual investigation has shown that smaller prominences tend to be more dy¬ 
namic than the larger ones, and often disappear from the limb much faster than 
what the solar rotation would imply. In addition, it was found that eruptive 
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(a) (b) 
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Figure 8. | (a) [ Distribution of the lengths of prominences. [(^Distribution of t he he ights above 
solar limb. |(c)| Distribution of the shapes in terms of length to height ratio. |(d)| Latitudinal 
distribution of prominences ; neg ative latitudes correspond to the s outh ern hemisphere, positive 
to the northern hemisphere |(e) j Distribution of centroid velocities. [(^Distribution of observed 
lifetimes. The parameters in panels a-e are averaged over the prominence lifetime. 


prominences show in general the same distribution of the inferred parameters, 
except for the fact that they tend to be among the bigger prominences. This 
suggests the existence of a certain critical mass, which should be reached before 
a prominence can become unstable. 

4.2. Estimation Models 

In contrast to determining spatial and dynamic characteristics, estimation of 
the mass and, consequently, the gravitational energy of prominences proves to 
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be a much more difficult task. The most precise methods developed so far are 
geuerally based ou measuriug the fractiou of backgrouud radiatiou absorbed 
by a promiueuce, most commouly iu the irou EUV hues (Gilbert, Holzer, aud 
MacQueeu, 2005 Schwartz et aL, 2012). Uufortuuately, this approach cau be 


ouly used iu a limited uumber of cases, wheu promiueuces are observed agaiust 
a sufficieutly uuiform backgrouud, which does uot chauge much over time. We 
have attempted, however, to apply this techuique usiug the 171 A chauuel data 
of TESIS, but as expected, with uo promisiug result. luterestiug developmeuts to 
this method have beeu proposed by Laudi aud Reale (2013) aud Williams, Baker, 
aud vau Driel-Gesztelyi (2013). These authors employed observatious iu several 
EUV hues with differeut absorptiou coefficieuts iu order to exclude reduudaut 
parameters aud thus to achieve higher precisious. Both these techuiques, how¬ 
ever, require promiueuce material to be co ol euo ugh so that it would uot emit iu 
the He ll 304 A hue. Low, Eoug, aud Eau (2003) have also proposed to estimate 
a promiueuce mass as the miuimal amouut of plasma ueeded to maiutaiu stable 
the curreut maguetic coufiguratiou; evideutly, this approach is uot applicable to 
our case either. Thus, we have uo better optiou thau to limit our study with ouly 
rough estimates of promiueuce masses aud euergies. Though beiug iuapplicable 
for the detailed studies of promiueuces, such estimates, however, cau be fairly 
useful for iuvestigatiug the geueral properties of a substautial uumber of eveuts. 

Iu our attempt to estimate the mass of a promiueuce, we rely ou the fact 
that the promiueuce plasma proves to be optically thick wheu observed iu the 
He II 304 A hue. This implies that the observed radiatiou is maiuly formed iu 
au outer layer of thickuess Lq, where the radiatiou from the uuderlyiug layers 
of the promiueuce is absorbed by helium ious. Erom the Beer-Lambert law, the 
thickuess Lq is related to the meau uumber deusity of partly iouized helium 
atouis di as 


Lq 1/e ^ l/riHeii , (4) 

where 5 is the meau absorptiou coefficieut of the emittiug layer aloug the hue of 
sight. We therefore assume that this outer layer is isothermal (with temperature 
To) aud homogeueous aloug the hue of sight, which meaus that the observed flux 
F is giveu by 

F = G{Tq)EM{Tq) ^ UHeii'^eLo ^ Uq , (5) 


where Ue is the electrou deusity, G is the coutributiou fuuctiou aud EM is the 
emissiou measure. Numerical simulatious by Gouttebroze aud Labrosse (2009), 
which cousider a much more realistic radiatiou trausfer, iucludiug scatteriug of 
iucideut radiatiou from the Suu, imply a more complex depeudeuce. However, 
for typical plasma temperatures aud deusities iu promiueuces these simula¬ 
tious suggest that a liuear depeudeuce as iu Equatiou (§ is au acceptable 
approximatiou. 

Iu the simplest case, oue cau relate the volume of a promiueuce, U, to its 
observed area, H, as V r\j A3/2. Here, we use au alteruative method iustead, aud 
estimate the thickuess of a promiueuce, T, at each pixel of the detected area 
assumiug that promiueuces have a more or less rouud sectiou. This meaus that 
the promiueuce attaius its maximum thickuess Dmax at poiuts that are the most 
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Figure 9. Total mass |(a)| and gravitational energy |(b)| of prominences simultaneously ob¬ 
served in a single image as they vary through the observation period. Grey dots correspond to 
individual images and the black line is the one-day smoothed trend. 


distant from the boundaries and that Dmax is equal to twice that distance. For 
the rest of the pixels, we calculate the thickness D as: 


D — 2 y {D^ax R) 5 


(6) 


where R is the minimal distance from the given point to the prominence edge. 
However, during our study, we have found that for the set of detected promi¬ 
nences this method gives a clear dependence V ^ ^i-33±o.07^ which is more 
consistent with the fact that prominences tend to be elongated rather than 
spherical in shape. Both methods, however, are far from being precise as they 
do not take into account projection effects and, again, can be used for a rough 
evaluation only. 

These two models allow us to obtain estimates of the mass and, therefore, of 
the gravitational energy of a detected prominence. We choose the proportionality 
coefficient in Equation <§ such that, for all the detected prominences, the mean 
value of the electron density would equal a typical value of 2 x 10^^ cm“^, which 
was found in previous studies of individual prominences (Bommier, Leroy, and 


Sahal-Brechot, 1986; Bommier et a/., 1994 Wiik, Heinzel, and Schmieder, 1994). 


We have to assume here, of course, that the electron density inside a prominence 
remains proportional to that in its outer shell along the line of sight. Earlier 
studies of prominences have shown that their ionisation degree is typically in the 
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Figure 10. Geometrical model used to evaluate the fraction of prominences visible off-limb. 
The prominence BT is located at a critical position, where it may become not visible above 
the limb. 


range of 0.5-0.8 (Orrall and Schmahl, 1980[ Kanno, Withbroe, and Noyes, 1981) 


and their helium abundance is around 0.1 (Heasley and Milkey, 1978). Therefore, 
suggesting that most of the prominence mass is due to the hydrogen atoms, and 
the low helium abundance will not imply a large error, we calculate the mean 
density of the prominence plasma by multiplying the inferred electron density by 
the proton mass. Subsequently, multiplying this plasma density by the previously 
calculated thickness of the prominence, we obtain the column mass at the every 
point of its visible surface. Thus, we can obtain the total mass of a prominence 
by integrating these column masses over the observed area of a prominence, 
and the gravitational energy by integrating the column masses multiplied by the 
gravitational potential (here, we assume that the detected prominences are in 
the plane of the sky). This is, by no means, a precise method, but with the data 
available, it still can be used for rough estimates of these important values. 


4.3. Mass and Gravitational Energy 

Having summed up masses and energies of all individual prominences detected 
in a single image, we can calculate the total mass of prominences observed 
simultaneously off-limb. During the observation period, the total mass of the 
prominences varied greatly with time between and kg (Figure [9^, 

as well as the total gravitational energy, which was found to be in the range 
from up to erg (Figure [91^. 

In effect, we can observe, and therefore detect, only a certain fraction of all 
prominences present simultaneously on the Sun, which are located close to the 
limb and are high enough to be seen off-limb at their position. To estimate this 
fraction, we develop a simple geometrical model presented in Figure Here, a 
prominence BT, which has a typical height ho, is located at a critical position 
where it becomes no more visible above the limb. We can assume ho — 40 Mm, 
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taking the most probable value from the distribution in Figure however, 
since prominences are not visible below ^15 Mm because of strong background 
radiation from the EUV spicules, we take here ho = 25 Mm and the solar radius 
Rq = 715 Mm. According to this model, we are able to observe only those 
prominences that are located between the planes APK and BQL (Figure 10). 
The area of the solar surface confined between these two planes is given by: 


Ayis = 2ttRq21 = 47ri?Q sin a , (7) 

where 21 = PQ and a is the NOB. Since ^ttRq is the total surface area of the 
Sun, we can observe off-limb a fraction of around 


sin a = — y^2ho/i70 ^ 0.25 (8) 

of all prominences present anywhere on the Sun, including those visible as fila¬ 
ments on the disk. Taking this into account, we conclude that at solar minimum, 
all quiet-Sun prominences contain simultaneously about lO^^-lO^^ kg of plasma 
and their total gravitational energy is of the order of lO^^-lO^^ erg. 

The big number of detected events enables us to investigate the form of their 
gravitational energy spectrum, which is shown in Figure |lla[ Although this 
spectrum may be shifted towards lower or higher energies due to the arbitrary 
choice of the coefficient in Equation its shape is only biased by our assump¬ 
tions about the prominence morphology and density distribution along the line of 
sight. The main feature of this spectrum is that it has a clear maximum at around 
1028.6 

erg. We believe that this behaviour can be explained by the existence of a 
low sensitivity threshold in our algorithm. Most likely, we are not able to detect 
small and low-lying prominences in the 304 A images, as they become barely 
distinguishable from the strong background below 15 Mm (Eignre[3a|). Indeed, we 
have detected no prominences lower than 20 Mm (EigurejS^. On the other hand, 
there can be no such small prominences at all, as they may become unstable at 
these spatial scales or indistinguishable from smaller objects seen above the limb. 
If so, the conventional definition of a prominence becomes no more applicable in 
this case. Anyway, we consider only that part of the spectrum above 10^^-^ erg to 
have a physical meaning. In a log-log plot (Eigure pM^ , it shows a clear power- 
law dependence, its index given by the slope of the distribution that we find 
to be —1.1 ± 0.2. Hence, in this case, the power-law index is close to —1, which 
means, that at least within the sensitivity range of our method, the gravitational 
energy is close to be evenly distributed among prominences of different sizes. 


5. Discussion 

As we mentioned above, most of the recent prominence studies are limited to 
particular cases, in general to the biggest events, which are easier to identify 
by eye. At the same time, the whole set of prominences is of great interest, 
and several questions remain, to which answers cannot be given by studies 
of individual events. It is unclear, for example, how the physical properties of 
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(a) (b) 



Figure 11. |(a)| Gravitational energy spectrum of the prominences. [(^ The slope of the same 
distribution for a limited energy range in a log-log scale. 


prominences vary from one to the other and what their long-term behaviour is, 
z.e., how they change throughout the solar cycle (for a discussion, see a review by 
Patsourakos and Vial, 2002). Detailed catalogues that give not only qualitative, 
but also quantitative information and comprise a significant number of their 
measurable parameters would enable us to perform large statistical studies of 
these phenomena. These studies, perhaps being not as precise as the individual 
ones, allow us, however, to cover a much wider range of events, maintaining a 
consistent point of view, less dependent on the selection of the objects examined. 

A catalogue of prominences above the limb would also be an appropriate com¬ 
plement to the existing catalogues of filaments on disk (Zharkova et a/., 2005), 
as they show prominences in two orthogonal projections. Equally promising are 
combined analysis of prominence catalogues with those of associated events in 
the solar corona. For example, deeper conclusions can be achieved by combining 


prominence catalogues and existing CME catalogues (Gopalswamy et a/., 2009 
Robbrecht, Berghmans, and Van der Linden, 2009 Floyd et a/., They are 


both observed in the direction orthogonal to the line of sight, and comparative 
studies of a great number of events would help to advance our understanding 
of the mechanisms underlying CME initiation. A joint study including flare 
catalogues (Watanabe, Masuda, and Segawa, 2012 Aschwanden et aL, 2014) 
would be of no less interest, as prominence eruptions are often an essential part 
these complex processes. 

Within the limits of this study, we have investigated a relatively short period. 
Nevertheless, it provides a lot of new information about general properties of 
the entire set of prominences, including rather small ones. One should also note 
that this observation period coincides with the end of an extraordinarily long 
solar minimum, which is a unique event by itself. The main advantage of the 
algorithm developed is a well-founded technique of background removal and a 
more accurate reconstruction of prominences, which increases the reliability of 
the measured prominence properties. Moreover, this algorithm is able to dis¬ 
criminate between the prominence and off-limb manifestations of active regions 
without the need of learning the method previously or using multi-wavelength 
observations. 
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The obtained values of mass and gravitational energy of prominences, despite 
the simplicity of the evaluation method used, are, nevertheless, the first attempt 
to determine these important properties for a substantial number of prominences. 
These data are useful to constrain existing models of prominence formation and 
evolution and to study their role in the general mechanisms of mass and energy 
transfer in the solar corona. It is also relevant that these properties are not only 
obtained for the largest events, but for the whole ensemble. Therefore, they give 
an early insight of how these parameters are distributed among prominences of 
different size. 

Because of the relatively short operation time of TESIS, our next goal is to 
adapt our algorithm to process the data obtained by other instruments, primarily 
those on board SDO. These data will permit us to examine prominences on much 
larger time scales and thus to detect long-term variations of their properties 
throughout the solar cycle, if any, and to find out why this solar minimum 
was so different from the others. For now, we can see two major problems that 
we will encounter in this endeavour. First, our study did not include active 
prominences, as they were not numerous during the examined period. Closer to 
the solar maximum, however, the number of this type of prominences will greatly 
increase along with the general increase of solar activity. Nevertheless, we believe 
that our algorithm can be successfully used to detect these prominences as well, 
against the enhanced active region background. The other problem is that our 
method of determining the background becomes inapplicable because of a more 
complex morphology of the corona during the rest of the solar cycle. Anyhow, 
the statistical approach used here remains a powerful tool to deal with this issue 
as well. 
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